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EXACT  DYNAMICALLY  CONVERGENT  CALCULATIONS  OF  THE 
FREQUENCY-DEPENDENT  DENSITY  RESPONSE  FUNCTION 

J .  Hong 

Department  of  Physics  Education,  Seoul  National  University,  Korea 

M.  Howard  Lee 

Physics  Department,  University  of  Georgia,  Athens,  GA  30602 


ABSTRACT 

A  general  expression  for  the  response  function  is  derived 
by  the  method  of  recurrence  relations.  Memory  effects  appear  as 
corrections  to  the  dynamic  RPA  form.  The  dynamic  structure  for 
the  3d  electron  gas  is  calculated  to  third  order  and  compared 
with  A1  data  at  a  large  wave  vector.  Also  shown  is  the  (dynamic  local- 
field  tern. 


PACS:  05.30  Fk 
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The  dynamic  response  in  strongly  interacting  homogeneous 


many-fermion  systems  has  drawn  considerable  attention  recently. 


Let  such  a  many-body  system  be  described  by  H  =  +  v,  where 


Hq  represents  the  kinetic  energy  and  V  the  interaction  energy. 


and  the  response  in  the  system  by  the  response  function  t 

where  k  and  ca  are,  respectively,  wave  vector  and  frequency  (V=l)  . 


It  is  well  known  that  the  response  function  may  be  put  in  the 


form 


Xk^*^)  =  /(I  +  (a>)) 


where  X]^(o)  (w'  is  the  response  function  due  to  and  Aj^(aj)  some 


unknown  function  of  V.  Various  dynamic  RPA  theories  are  equivalent 


to  taking  Aj^(w)=Aj^(0)  iAj^.  For  example,  Aj^  =  Vj^  gives  the  simple 

RPA  where  Vj^  is  the  two-body  interaction,  and  Aj^  “  ”  ®k^ 

gives  the  generalized  RPA  where  Gj^  is  a  local  field  term  effective 

3  4 

over  the  Thomas-Fermi  length.  '  There  are  large  numbers  of  RPA 


based  studies  and  they  have  been  routinely  used  to  interpret 


experiments . 


The  dynamic  RPA's,  however,  can  break  down  especially 


in  the  regime  of  metallic  densities  at  large  wave  vectors.  Clearly, 


one  needs  to  restore  the  frequency  dependence  e.g.,  Aj^(a))=  Vj^  (1  -  Gj^(( 

_ ^ _ _  _ £ _ ^  ^  _ f _ ^  6—9 


Several  people  have  obtained  asymptotic  conditions  for  • 


Otherwise,  there  are  at  present  no  exact  general  expressions  known. 


For  a  3d  Coulomb  gas  Devreese  et  al^^  and  Holas  et  al^^  have  calcu¬ 


lated  by  some  approximate  techniques.  Unfortunately,  they 

find  in  their  calculations  several  unphysical  divergences.  This 


sort  of  situation  has  limited  our  understanding  of  the  dynamic 


structure  of  metals  at  large  wave  vectors.  In  particular,  whether 


•-'I*.;* 


the  observed  fine  structure  arises  from  multi-pair  excitations 

4 

or  from  some  other  strong  correlations  remain  unresolved. 

In  this  Letter,  we  obtain  a  general  expression  for  Aj^(u)  by  the 
method  of  recurrence  reltions  and  perform  exact  dynamically- 
convergent  calculations  based  on  it.  Our  results  are  compared 


with  A1  data. 

The  time  evolution  of  the  density  f^ctuation  operator 
may  be  given  an  orthogonal  expansion  viz., 

where  |  |  forms  a  complete  set  of  basis  vectors  spanning  the 

d-dimensional  Hilbert  space  of  and  are  autocorrelation 

functions.  There  exist  recurrence  relations  RRI  and  RRII, 


respectively,  for  f^'s  and  a^*s: 


■v+1 


+  K  *v-l 


(2) 


(3) 


where  O  <  v  ^  d-1,  f^=  ilH,  f^] ,  a^=  da^/dt,  =  (f^r f^) / (f^_l » 
where  the  inner  product  means  the  Kubo  scalar  product  and  by 
definition  f_2^=0,  a_^  =  O  and  A^  =  1.  According  to  the  method 


of  recurrence  relations ,  one  obtains  the  basis  vectors  f ^ ' a  by 


RRI  (2)  from  which  the  recurrants  A^*s.  One  can  then  deduce  the 


autocorrelation  functions  a  's  by  RRII  (3)  and  obtain  the  time 
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evolution  of  Pj^. 

If  we  choose  f^  =  Pj^,  then  by  linear  response  theory 


ao(t)  »  (pj^{t),  '  which  is  the  relaxation  function, 


and  A^aj^(t)  Xj^  = 

Aiai(t)  *  -a^(t).  (4a 

By  applying  the  Laplace  transform  operator  T  and  with  a^(t»0)  * 


3 


we  obtain 


Ajj^a^(z)  *  1  -  z  a^(2) 


(4b) 


where  a^(z)  =  T(a^(t)l.  Hence,  from  (4b), 


a^(z)  =  (z  +  A^b^(z)) 


where  b^(z)  =  aj^(2)/a^(z) ,  Now  combining  (4b)  and  (5)  and  using 
the  identity  Aj^aj^(z)  =  x(2)/x»  suppressing  k-dependence  altogether. 


we  get 


Aj^b^(z) 
z  +  A*bj^(z) 


Since  eq.  (7)  is  valid  for  amy  Hermitian  model,  it  also 
applies  to  (ideal  version).  Now  consider  * 

Since  =*  f^  from  (2),  we  have  =  i[H,Pj^]  = 


Hence , 


J: -  =  X,-.  , 

We  divide  (6)  by  its  ideal  version  and  obtain 


.X  (?) 

X  ^°\z) 


*1  ^  (z)  r  *  Aj^bj^(z) 
(2)  L  2  +  Ai®^bJ 


(z) 


1  i  A(z)x''"'  (2) 


(8a) 


(8b) 


Observe  that  (8b)  is  in  the  form  of  (1) .  By  using  the  ideal 
version  of  (6) ,  we  find  from  (8a,b)  the  desired  expression: 

A(z)  =  I(x)"^  -  [(bj^(2))~^  -  (b^°^(z))"^l 

=  A  +  X(z).  (9) 

The  first  bracketed  term  A  is  z- independent.  Thus,  the  BPA 

theories  are  valid,  i.e.,  A(z)  =  A,  if  and  only  if  b^^  (z)  =  b^  (2 

We  shall  now  examine  thr  z-dependent  part.  EarlierT 

b^(z)  was  introduced  in  place  of  a^(z)/a^(z).  According  to  the 
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method  of  recurrence  relations  there  is  actually  a  whole  family 

of  b^(t)  =  T  ^[b^(z)],  1^  V  <d-l,which  define  the  time  evolution 

d-1 

of  the  generalized  random  force  F.  for  p.  as  P.  (t)  =  J  • 

K  K  K  V  I  ^ 

The  random  force  lies  in  a  linear  manifold  of  the  Hilbert  space 
of  pj^.  The  autocorrelation  functions  b^'s,  sometimes  referred 
to  as  the  memory  functions,  also  satisfy  a  recurrence  relation: 


+  Vi<‘> 


<10) 


where  b^(t)  =  0  and  l<v<d-l.  That  is,  for  example, 

bi(t)  =  bj^(A2,A3, - t)  .  (11) 

,  »  c/ne  fc  tnhr<ic  •i-tiitr' 

Hence,  if  \  2  jc  y»  ^  d-1,  memory  effects^are  always  present 

in  the  response  function  and  are  manifested  through  z  dependence 

inA(z).^^ 


The  recurrants  A^'s,  which  are  relative  norms  of  the  basis 
vectors,  are  model-dependent.  They  are  basic  elements  of  the 
dynamic  structure  and  their  form  shapes  the  time  evolution.  In 


some  special  physical  regimes  of  certain  memy-body  models  they 
can  be  calculated  to  any  order  and  one  can  use  them  to  obtain 
the  memory  function  from  RRII  (10) .  For  noninteracting  systems 
e.g.,  an  ideal  3d  electron  gas,  the  recurrants  are  in  effect 
kno%m  to  all  orders  at  T  =  0  since  (z=iw)  is  analytically 
given  in  the  ground  state.  But  for  interacting  systems  generally 
only  the  first  few  orders  of  the  recurrants  have  been  accurately 
calculated  presently.  Hence,  one  C5mnot  obtain  their  memory 
function  by  (10)  and  one  cannot  use  the  general  expression  (9) 
to  explicitly  calculate  e.g. ,  im  metal  densities  at 

large  wave  vectors.  We  propose  here  a  technique  by  which  one 
C2un  use  (9)  to  calculate  A  (z)  term-by'>term  up  to  the  available 
order  of  the  recurrants  given  that  the  "ideal”  recurremts  are 
known  to  all  orders. 

From  (9)  and  (11)  we  see  that 

X  =  ({a^};{a^^°^};z),  v  >  2  (12) 

assuiaing  d  -  «>.  Hence  we  can  write 
X  =  lim  X 

n 

n  •►  «»  (13) 


where 


Xn  =  X(A2...Ajj 


(14) 


The  process  of  replacing  X  by  X  ^  has  the  following  significanceT 
If  X  is  used  in  place  of  X  in  (9),  the  resulting  x(z)  satisfies 
the  frequency  moment  sum  rules  exactly  up  to  and  including  the 
2n+l  st.  It  satisfies  the  higher  moments  to  the  accuracy  of 
the  substitution  of  A  by  m>n+l.^® 


6 


By  this  process  one  can  obtain  an  expression  for  the  response 
function  in  terms  of  the  recurrants  which  is  exact  to  a  given 
order.  One  can  continue  this  process  to  the  highest  available 
recurrant.  An  infinite-order  expression  evidently  is  equivalent 
to  an  exact  expression  for  the  response  function. 

Using  (9)  euid  RRII#  we  can  systematically  obtain  e.g. 

0  (15) 

X^^^X3  =  (n2  -  n3  R(2))Q(z)/(1  +  113  R(z)) 

etc.  where  rij  =  j  =  2,3 

Q(z)  =  . 

R(z)  =  1  -  Y  Q(z),  Y  =  ^ 

Observe  that  X  X^^  gives  the  generalized  RPA.  For  k  z  kp  where 
kp  is  the  Fermi  vector,  multipair  excitations  exist  which  are 
contained  in  A2A 3  etc.  They  can,  therefore,  contribute  to  the 
response  function  via  frequency  dependence  beyond  first  order. 

The  convergence  of  our  term-by-term  calculations  should  be 
relatively  rapid  since  our  procedure  uses  infinite  continued 
fractions  at  all  stages  which  are  known  to  give  good  convergence.^^ 
Finite  order  calculations  are  meaningful  if  x^°^  0. 

To  illustrate  our  technique  ,  we  calculate  the  dyneunic 
structure  factor  Sj^(u)=-7t”^  Ira  x(z“i“)  the  3d  electron  gas 
to  third  order,  the  highest  order  possible  based  on  the  presently 
available  recurrants  which  are  A^A2A3.  For  our  calculations  we 


k  —  1.6  kp  and  —  2.0  (c£.  A1  ®  2.07)  .  For  'these  values 
we  find:  =  0.1297,  n3  =  0.0364,  and  y  =  0.5619. 


Observe  that  in  our  finite  order  calculations  Sj^(a))=0 


for  tii  >  5.76  £_. 
-  F 


Shown  in  Fig.  lA  is  the  first-order  corresponding  to  the 
generalized  RPA  (dotted  line) .  It  shows  almost  no  fine  structure. 
The  second  order  (dash-dot  line)  and  the  third-order  (solid  line) 
begin  to  show  some  structure.  The  calculated  amplitudes  are  all 
absolute , 

1  O 

In  Fig.  IB  the  third-order  dynamic  structure 

(solid  line)  is  compared  with  the  simple  RPA  (dotted  line)  and 

a1  data  (dashed  line) .  The  experimental  amplitude  is  adjusted 

to  coincide  with  our  third-order  calculated  one.  The  simple 

RPA-calculated  amplitude  is  absolute.  We  observe  that  the  third- 

order  calculated  structure  factor  shows  some  resemblance  to  the 

experimental  curve  especially  in  the  shoulder.  The  peak  position 

is  much  closer  to  the  experimental  peak  position  than  that  of 

20 

the  simple  RPA. 

Using  the  definition  ^(w)=  Vj^  (l-Gj^(a))  )one  can  also  extract 


the  corresponding  frequency-dependent  local  field  term  Gj^(b)). 

These  results  are  shown  in  Fig.  2.  The  first-order  (dotted  line) 
is  absent  in  Im  Gj^((u)  since  1^^  =  0.  The  second-order  (dash-dot 
lines)  and  third  order  (solid  lines)  are  well-behaved,  containing 

no  infinities.  In  Re  G,(w)  there  are  cusp-like  peaks  at  u  -e— 0.-5- 

k  .  10  / 

and  It  is  interesting  to  note  that  Devreese  et  al  ^  Hclca  ci 

encountered  divergences  in  their  calculation  of  Gj^  (w)  at  ^tpproxiraately 
these  frequencies. 


This  work  was  supported  in  part  by  DOE,  ONR,  the  Korea  Science  and 
Engineering  Foundation,  and  the  Research  Institute  for  Basic  Science 
Seoul  National  University. 
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To  obtaim  the  high  frequency  tail  (see  Fig.  IB)  one  must 
consider  an  infinite  order  calculation. 
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Figure  Captions 

1.  Dynamic  structure  vs  frequency.  A.  first-order  (dotted), 
second-order  (dash-dot),  third-order  (solid).  B.  experimental 
(dash) ,  simple  RPA  (dotted) ,  third-order  (solid) .  E^Fermi 
energy 

2.  Dynamic  local  field  correction  vs  frequency.  First-order 
(dotted) ,  second-order  (dash-dot) ,  third-order  (solid) . 

Ej,  Fermi  energy 


